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Abstract 

The four exactly-solvable models related to non-sinusoidal coordinates, namely, the Coulomb, 
Eckart, Rosen-Morse type I and II models are normally being treated separately, despite the 
similarity of the functional forms of the potentials, their eigenvalues and eigenfunctions. Based on 
an extension of the prepotential approach to exactly and quasi-exactly solvable models proposed 
previously, we show how these models can be derived and solved in a simple and unified way. 
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I. INTRODUCTION 



Exactly-solvable models are important in any branch of physics. They allow a complete 
understanding of the dynamics of the corresponding systems, and thus serve as paradigmatic 
examples for the field of studies concerned. However, exactly-solvable systems are rather 
scanty, and hence any means to find them are always welcome. It is thus very interesting that 
most exactly-solvable one-dimensional quantum systems can be obtained in the framework 
of supersymmetric quantum mechanics, based only on the requirement of shape invariance 
[1]. Recently, the number of physical systems which we can study analytically has been 
greatly enlarged by the discovery of the so-called quasi-exactly solvable models [2-7] . These 
are spectral problems for which it is possible to determine analytically a part of the spectrum 
but not the whole spectrum. 

In [8] a unified approach to both the exactly and quasi-exactly solvable systems is pre- 
sented. This is a simple constructive approach, based on the so-called prepotential [9, 10], 
which gives the potential as well as the eigenfunctions and eigenvalues simultaneously. The 
novel feature of the approach is the realization that both exact and quasi-exact solvabilities 
can be solely classified by two integers, the degrees of two polynomials which determine 
the change of variable and the zero-th order prepotential. Most of the well-known exactly 
and quasi-exactly solvable models, and many new quasi-exactly solvable ones, can be gener- 
ated by appropriately choosing the two polynomials. This approach can be easily extended 
to the constructions of exactly and quasi-exactly solvable Dirac, Pauli, and Fokker-Planck 
equations. 

The exactly solvable models that are generated by the prepotential approach in [8] 
are related to the so-called sinusoidal coordinates [11]. These are coordinates z{x) whose 
derivative-squared z''^{x) (henceforth the prime denotes derivative w.r.t. the variable x) is 
at most a quadratic polynomial of z (or equivalently, z" is at most linear in z) . Such coordi- 
nates include quadratic polynomial, exponential, trigonometric, and hyperbolic types. The 
corresponding exactly solvable models are the six systems hsted in [1]: the shifted oscillator, 
three-dimensional oscillator, the Morse, Scarf type I and II, and generahzed Poschl-Teller 
potentials. However, the other four models in [1], namely, the Coulomb, Eckart, Rosen- 
Morse type I and II models cannot be generated by the approach as presented in [8] . These 
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four models are based on a change of coordinates which are non-sinusoidal,^ which accord- 
ing to the discussions in [8] can only be quasi-exactly solvable. This is consistent with the 
Lie-algebraic approach to quasi-exactly models, as these four models, unlike the other six 
systems, cannot be obtained as the exactly solvable limits of some quasi-exactly solvable 
systems [3, 4]. 

The potentials of these later four systems are usually given as follows [1]: 

Coulomb : A (A - 1) ^ - 25^; 

Eckart : A(A ^ a) cosech^ ax — 2i?coth ax; 

(1) 

Rosen-Morse II : —A {A + a) sech^ ax + 2Stanh ax; 
Rosen- Morse I : A{A — a) cosec^ ax + 2Scot ax, 

where A, B and a arc real constants. Despite the similarity of the functional forms of the 
potentials, their eigenvalues and eigenf unctions, these models are normally being treated 
separately in the literature. 

The purpose of this paper is to show how these four systems can be generated by a 
simple extension of the prepotential presented in [8]. What is more, they can be treated in 
a unified manner. Put simply, the required extension is simply to allow the coefficients in 
the prepotential, which are assumed to be constants in [8], to be dependent on the number 
of nodes of the wave functions in such a way that the coefficients of all terms involving z in 
the potential are real constants. 

To bring out the close similarity of these four potentials, we find it convenient to express 
the functions cosech, sech and cosec in terms of coth, tanh and cot, respectively, in (1). The 
resulted forms of the potentials, the relevant change of coordinates, and their eigenvalues 
are listed in Table 1. For clarity of presentation, we adopt the unit system in which h and 
the mass m of the particle are such that h = 2m = 1. Also, without loss of generality, we 
have absorbed the scale factor a into x, or equivalently, we set a = 1. 

Prom the table, we see that these four models involve a change of variable z{x) whose 
derivative is of the form 

/ = A - (2) 

^ We note here that the Coulomb potential can in fact be treated using sinusoidal coordinate, which we 
shall present in the Appendix. 
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TABLE I: The four exactly solvable models based on non-sinusoidal coordinates. The potential 
V{x), its relevant non-sinusoidal coordinate z{x), the derivative z'{x), and the eigenvalues Ej\i 
{N = 0,1, . . .) are listed. Without loss of generality, we absorb a into x, or equivalently, we set 
a = 1. The range of a; is: x E [0, oo) for the Coulomb and Eckart potentials, x G (— oo, oo) for the 
Rosen- Morse II potential, and x G [0, tt] for the Rosen-Morse I potential. 



V{x) 



z{x) z'{x) 



En 



Coulomb 



A{A-l)^-2Bl 



—z 



Eckart A{A - 1) coih'^x - 2Bcoth x coth x I - z^ -j^^-A{2N + \)-N'^ 



Rosen-Morse II A ( A + 1) tanh^x + 2Btanh x tanh x \-z^ - + A {2N + 1) - iV' 



Rosen- Morse I A{A-1) cot^x -|- 2Bcot x cot 



X 



-1-z' 



B2 



(A+N)' 



+ A{2N + 1) + iV2 



0, 
1, 



-1, 



Coulomb; 

Eckart, Rosen-Morse II; 
Rosen- Morse I, 



and that the potentials can be cast in the same form 



(3) 



V{x) = A{A - l)z\x) - 2Bz{x). 



(4) 



In this form, the four potentials are regular in the variable z, with singularity only at z — oo 
and/or z = — oo. As the function z'^ is a forth-degree polynomial in z, the coordinate z{x) is 
non-sinusoidal. Naively any system based on such coordinate will be quasi-exactly solvable 
according to the discussion in [8]. In the rest of the paper, we shall show how the exactly 
solvable models in (4) can be obtained from (2). 
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The organization of the paper is as follows. In Sect. II we review the main points of 
the prepotential approach of [8] relevant to our present discussions. Sect. Ill shows how 
the exactly solvable models (4) can be derived from the non-sinusoidal coordinates (2) by 
a simple extension of the approach in [8]. The four specific models are then discussed in 
detail in Sect. IV to VII. Sect. VIII concludes the paper. In the Appendix, we show how 
the Coulomb system is treated based on sinusoidal coordinate. 



II. PREPOTENTIAL APPROACH 

The essence of the prepotential approach is as follows. Consider a wave function (p^ix) 
{N: non-negative integer) which is defined as 

0^(x) = e-^°(-W(^), (5) 

with 

' 1, 7V = 0; 



Pn{z) = < 



(6) 

ULii^-Zk), N>0. 



Here z = z{x) is some real function of the basic variable x, Wo{x) is a regular function of 
z{x), and z^s arc the roots of Pn{z)- The variable x is defined on the full line, half-line, or 
finite interval, as dictated by the choice of z{x), We have assumed that the only singularities 
of the system are z — oo and/or z — — oo, as is the case for the four models of concerned 
here. The function Pn{z) is a polynomial in an {N + 1) -dimensional Hilbert space with the 
basis (1, z^, . . . , z^). Wo{x) defines the ground state wave function. 
We rewrite ^at as 

0JV = exp {-Wn{x, {zk})) , (7) 

with Wn given by 

N 

Wm{x, {zk}) = Wo{x) -J2ln \z{x) - Zk\. (8) 

k=l 

Operating on by the operator —(P/dx^ results in a Schrodinger equation Hi>j(t)j^ — 0, 
where 

Hn^-^,+V^, (9) 

= w'^ - w;,. (10) 



Hence the potential Vn is defined by Wn, and we shall call Wn the A'"th order prepotential. 
Prom Eq. (8), one finds that Vjv has the form Vn — Vo + AVn- 



Thus the form of Vn, and consequently its solvability, are determined by the choice 
of Wo{x) and z'"^ (or equivalently by z" — {dz'"^ / dz) / 2) . Let Wqz' and z''^ be taken as 
polynomials in z. In [8], it was shown that if the degree of W^z' is no higher than one, and 
the degree of z'"^ no higher than two, then in Vn{x) the parameter N and the roots ZkS, 
which satisfy the so-called Bethe ansatz equations (BAE) to make the potential analytic, 
will only appear in an additive constant and not in any term involving powers of z. Such 
system is then exactly solvable. If the degree of one of the two polynomials exceeds the 
corresponding upper limit, the resulted system is quasi-exactly solvable. 

Now we are interested in constructing a system based on a transformed variable z{x) 
which is a solution oi z' — \ — z^ . Since in this case the degree of z'"^ is four, a direct 
application of the arguments in [8] would mean that no exactly-solvable system can be 
generated with whatever choice of W^q- 

But it turns out that with a slight extension of the methods in [8], one can generate the 
four potentials in the prepotential approach. The main observation is this. The classification 
of solvability given in [8] is vahd as long as the coefficients of the powers of z in W^z' are 
constants which are independent of N . If one allows these coefficients to be A'"-dependent, 
then it may be possible that in the Vn obtained, all the coefficients of terms involving powers 
of z are A^-independent constants, thus giving rise to an exactly-solvable system. This is 
in fact the case for the four systems mentioned above. We shall demonstrate this in what 
follows. 

III. EXACTLY-SOLVABLE MODELS WITH NON-SINUSOIDAL COORDI- 
NATES 

Let us choose z from a solution oi z' — \ — z^, and take 



Til 

1 




(11) 



W'^ = A^z + Ao, 



(12) 



6 



where Ai and Aq are real parameters. With this choice of Wq and 2;'^, we obtain from (11) 



Vo = ^1(^1 + 1)^2 + 2^1^0-2 + Al- XAi 



and 



AVn = -(^X-z^)< 
Using the identities 



fc=i ^ fc=i ^ 



E 

=1 

1 



y — 

AT 



N 



2E 



1 



E 



z 



k^l 

N 1 / „2 



2E 



k^l 



+ N{N-1), 



we rewrite AV^ as 

I ■'^ 1 

AVn ^-{X-z'')In (2A^ + n+1) + 2^ — 



N 

k=l ^ ~ 



{Ai + 1) Zfe + Ao + E 



l^k 



To remove the poles in AVn, we must demand that z^s satisfy the BAE 
^^i^ + (Ai + l)^, + Ao = 0, A; = l,2,...,iV. 

With these 2;fe's, only the first term in AVn remains, and the potential Vn becomes 



(13) 



• (14) 



(15) 
(16) 

(17) 



(18) 



(19) 



Vn = {Ai + N){Ai + N + 1) z^ + 2AiAoz + AI-X [{2N +l)Ai + N{N+ 1)] . (20) 

Now it is seen that N appears in the coefficient of z'^ term, and thus Vn represents a quasi- 
exactly solvable system, if Aq and Ai are some fixed constants. 

But in this case it is easy to obtain an exactly solvable potential. Let us choose Ai and 
Aq such that the combinations A = —{Ai + N) and B = —AiAq are A^-independent real 
constants. The choice of the signs in the definitions of A and B are for convenience and are 
not essential for the moment (they will have to be determined by the normalizability of the 
wave functions for physical systems). Consequently, Ai and Aq depend on N: 

B 



A,^-(A + N), Ao 



A + N' 



(21) 



Then the potential Vn becomes Vn{x) — V{x) — En, where 

V{x) = A (A - 1) z\x) - 2Bz{x), 



(22) 



as advertised in Sect. I, and 

B 



En — o — A 



A {2N +1) + N^ . (23) 



Now V{x) is independent of N, and can be taken to be the potential of an exactly-solvable 
system, with eigenvalues En {N = 0, 1,2, . . .). The corresponding wave functions (pN are 
given by (5) together with (12) and (21): 

0Ar~e(^+^)r'^^"(^)-^>A,(a;), iV = 0, 1, . . . (24) 

Prom (19) and (21), the BAE satisfied by the roots ZkS are 

J2^^-{A + N-l)z, + ^^^0, k^l,2,...,N. (25) 

In the sections below, we will show how the four potentials mentioned emerge from (22) 
with different A. But first a comment on the choice of the form of V{x) is in order. One 
notes that there is an ambiguity in the definitions of V{x) and En in Vn = V — En'- Vn is 
invariant under V ^ V — a and En — ^ En — a for any real a. This amounts to the choice 
of the zero point of V{x). The form V^susy ^(Jopted in supersymmetric quantum mechanics 
(e.g., in [1]), corresponds to the choice a = £'o so that the ground state energy is zero, i.e. 
£■0 = 0. In our case, y^usY -g obtained from the zero-th order prepotential Wo{N — 0) with 
N — (remember now that Aq and Ai in Wq depend on N): 

W;,{N = 0) = -Az + |, 

ySUSY ^ |^^2(^ ^ _ y^.^j^ ^ (26) 

= A{A - l)z^ -2BZ + — + XA. 



The energies are 



A^ 



B^ B^ 



A^-(A + NY + — ^. (27) 

^ ^ ^ A'' (A + Nf ^ ^ 

In the literature on supersymmetric quantum mechanics, Wq{N = 0) is generally called the 
superpotential. 
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IV. A = 0: THE COULOMB CASE 



In this case, A = and z' — —z^. A solution Ss, z — Xjx. The domain of x is x e [0, oo). 
The potential is 

n.) - ^ - 2^. (28) 

X 

From (24) the wave function is 

<^n{x) - x^+^e-^^ vn{z). (29) 



For (/)jv— >Oata; = and oo, one must have A > and S > 0. The eigen-energies are given 
by (23) 

En^-^£^,^ at. 0,1,... (30) 

The roots ^r^'s satisfy the BAE (25) with A = 0. It should be emphasized that the above 
expressions are valid for any real A > (not necessary integer) and i? > 0. 

The ordinary Coulomb potential is obtained by setting A = / + 1 (/ = 0, 1, . . .) and 
B — e^/2, where I is the orbital angular quantum number and e the electric charge. This 
gives the well known form of the Coulomb energies — — e^/4(A" + I + 1)^. To write the 
wave functions (29) in the familiar form, we express z and Z]^ in Pn{z) in terms of a; = 1/z 
and Xk = l/zk- This gives (by factoring out all z's and z^s) 

2 N 

(t)N{x) ~ x'+^e~2(ivW)^ n - ^k) , AT = 1,2,... (31) 

k=l 

Recall that pn{x) — 1 ior N — 0. The BAE satisfied by x^s are obtained from (25) by 
setting A = and changing Zk,zi to l/xk,l/xi. The final result can be written as 

y^^ + ^-^ = - A; = 1,2,..., AT. (32) 

^^Xk-xi Xk 2{N + l + iy ' ' ' ^ ^ 

Letting y = e^xj {N + / + 1), we can simplify the BAE to 

T.^— + — = i k^l,2,...,N, (33) 

where 7 = 2(i + 1). Eq. (33) is just the set of equations satisfied by the roots of the Laguerre 
polynomials L]^^(y) [9, 12]. Hence in terms of the variable y, the wave functions are 

My) - Z/'+^e-U^+i ^ (34) 
which are exactly the form given in Table 4.1 of [1]. 
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V. A = 1: THE ECKART CASE 



We take z — coth x. The range of x is again the half-hne x e [0, oo) . The potential, 
energies and wave functions are 

V{x) = ^ - 1) coth^a; - 2Bcoth x, 
En = o -Ai2N + l)- N^, (35) 

4>N ~ (sinh a;)^"^^e"^^ Pat (^) . 

Now the boundary conditions 0Ar— >Oasx— >0 and x — > oo require A > and B > {A+N^, 
respectively. Hence B must be greater than i.e. B > A^, in order to admit at least one 
bound state (corresponding to A?^ = 0). For fixed A and B > A'^, the maximal value of N is 
such that B > [A + N)"^ remains valid. Thus the number of bound states is iV + 1. 

Since zl ^ 1 one can divide (25) (with A = 1) by — 1. Writing the result in partial 
fractions, we obtain 

J^^^kl^^imil^O, *^l,2,...,iV, (36) 

- Zl Zk - 1 Zk + i. 

with 

a^-A-N+-^^, j3 = -A-N (37) 
A+N A+N ^ ' 

One recognizes that (36) are the equations satisfied by the roots of the Jacobi polynomial 
Pj^'^\z) [9, 12]. Hence the wave functions for the Eckart potential are 

0^ - (sinh a;)^+'^e-^^ p!^''^\z). (38) 
It is easy to check that (f)N can be expressed as 

M^) -^{z-l)Hz + \)- Pt^\z). (39) 
This is the form presented in [1]. 

VI. A = 1: THE ROSEN-MORSE II CASE 

Here A = 1 is the same as for the Eckart potential, but now we take a different solution 
z{x) = tanh x. This choice implies the variable x is defined on the whole line, — oo < a; < oo. 
The wave functions are 

(f)N ~ (cosh x)^+^ e-^^ pn (z) . (40) 
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The boundary conditions 0jv— >Oas|a;|— >oo lead to 



A + N<0, \B\<{A + Nf. (41) 

So in this case A must be negative, and B can have any sign as long as the inequality in 
(41) for B is satisfied. For easy comparison with the corresponding expressions in [1], we 
change A — > —A and B — > —B. With these new definitions of A and B, the potential (22) 
and energies (23) are written as 

V{x) = A{A+1) tanh^x + 2manh x 
En = -JX^^ + A{2N + 1)- N\ (42) 
^>0, \B\<{A- Nf . 

For fixed A and B the maximal value of N is such that \B\ < {A — iV)^ remains valid. The 
roots ZkS now satisfy a BAE which is the same as that in the Eckart case, but with the 
changes A — > —A and B ^ —B. Hence the wave functions are 

(t>N ~ (cosh x)-^^-""^ e"^^ Pp^\z), (43) 

where 

B B 

a^A-N+- , d = A-N-- . (44) 

A-N' ^ A-N ^ ' 

Again, one can rewrite ^at in the form given in [1] 

M^) ~ (1 - z)^ (1 + z)^ P'^'^\z). (45) 

We note there that in this case the two models defined by only a difference in the sign 
of B are simply mirror images of each other: they are related by the parity transformation 
X — > —X. 

VII. A = -1: THE ROSEN-MORSE I CASE 

Finally we come to the Rosen- Morse potential defined with the solution z = cotx for 
A = — 1. The system is defined on the finite interval x e [0, tt]. From the wave function 

0JV - (sina;)'^+'^ e"^^ pn (z) , (46) 
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one infers that A > while B is arbitrary. Again in order to facihtate comparison with the 
expressions in [1], we change B ^ —B. The two systems differing only in the signs of B are 
equivalent by reflection and periodicity. 
The potential and energies read 



V{x) = l)cot^x + 2Scotx 

B'' 



d2 

En = —2 +A{2N+l) + N\ (47) 



^^^-(A + N-l)z,--^^0, k^l,2,...,N. (48) 



The BAE (25) is 

i^k^k-zi ^ ' A + N 

This set of BAE is related to the Jacobi polynomials, but with imaginary argument. To see 
this, let us rewrite (48) in terms oly = iz. One gets 

5:^,i(a±i),i(^.0, k^l,2,...,N, (49) 

where 

a^-A-N-i-^. p^-A-N + i-^. (50) 
A+N A+N ^ ' 

Hence the wave functions are 

4>N - (sin xf^"" Pi^'^^ (iz) . (51) 

This is consistent with the expression in [1] in the form 

<Pn ~ (z' + l)"^ e^^'^*"'^ P^^^ {zz) . (52) 



VIII. SUMMARY 

In [8] a new approach to both exact and quasi-exact solvabilities was proposed. In this 
approach the solvability of a onc-dimcnsional quantum system can be solely classified by two 
integers, the degrees of two polynomials which determine the change of variable and the zero- 
th order prepotential. It was shown that exactly solvable models can only involve a change 
of variable which is sinusoidal, otherwise the system is quasi-exactly solvable. As such, 
four out of the ten exactly solvable models classified by shape invariance in supersymmetric 
quantum mechanics [1], namely, the Coulomb, Eckart, Rosen- Morse type I and II models, 
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are not covered by the approach as presented in [8]. In this paper, we have shown how 
these four models could be easily generated in a unified way in the prepotential approach 
by allowing the coefficients in the prepotential, which are assumed to be constants in [8] , to 
be dependent on the number of nodes of the wave functions. 

Thus with the results presented in this paper and those in [8], all the well known one- 
dimensional exactly solvable Schrodinger models have been generated by the prepotential 
approach in a way which, we believe, is much simpler and direct than other approaches. 
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APPENDIX: COULOMB POTENTIAL IN SINUSOIDAL COORDINATES 

Of the four potentials discussed in the main text. Coulomb potential can be treated in 
terms of sinusoidal coordinate, namely, z{x) = a; (as in the main text, for simplicity, we 
absorb any scale factor into x). Unlike those discussed in the main text, in the present 
case what we want is to generate a potential that is singular at x — 0. According to [8], 
the prepotential Wo{x) should be modified to Wo{x) — Wo{x) — A\-a\x\, where A is some 
real parameter, to ensure proper behavior of the wave function at the singular point x = 0. 
Naively, from the discussions n [8] , the degree of Wqz' should not be more than one in order 
to get exactly solvable model. So let us try Wq{x) — b {z' — 1). Thus the ground state wave 
function is 0o ~ exp(— Wq) ~ x^ exp{—bx). Normahzability of 0o requires A> and b > 0. 

Replacing Wo in (11) by Wq, we obtain 

X^ X 



13 



and 



N 



N 



^ k=i X Xk yi^^ 

Simple poles are removed if x^'s satisfied the BAE 



+ ^-6 



Xk Xi Xk 



(A.2) 



E 

l^k 



^ + — -6 = 0, k = l,2...,N. 



Xk Xi Xk 



Summing over k in (A. 3) gives 



AT 

A E bN^O, 

k=i 



and hence AV/v = —IbN/x. Finally, we have 

A{A-l) UA + N) 



Vi 



N 



- 2- 



+ 6^ 



(A.3) 



(A.4) 



(A.5) 



x^ X 

Since N appears in the 1/x term in V/v, the system is quasi-exactly solvable in this form. 
But as discussed in the main text, one can make this system exactly solvable by allowing b 
to depend on N in such a way that the coefficient of 1/x term is A^-independent, i.e. 

B 



A + N' 

where S is a real constant. So the Coulomb potential is given by 

A{A - 1) B B^ 



(A.6) 



2- + 



x^ 



\2 ■ 



X {A + NY 

This is consistent with the results given in (28) and (30). The BAE (A.3) becomes 

1 A B 



(A.7) 



E 



^A_ 



:, A; = 1,2...,A^. 



(A.8) 



i^k ^k-xi Xk A + N 

Letting A — I -\- 1 and B — e^/2, we recover the BAE in (32), and subsequently the wave 
function 0jv in (34). 
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